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Abstract 

The effective action describing the long range fluctuations in the high temperature phase 
of the electroweak standard theory is a strongly coupled SU(2)-Higgs-model in three di- 
mensions. We outline in detail a model in which the spatial correlation scales in this phase 
are calculated as inverse relativistic bound state masses. Selection rules for these states 
are derived. The correlation masses are calculated by evaluating the bound state Green's 
function. The scalar-scalar-potential and its influence on the masses is investigated. The 
predictions for the correlation masses agree very well with the lattice data available now. 
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1 Introduction 



The perturbative treatment of the electroweak standard model (SM) leads to a spectacu- 
larly good agreement with experiments. It works because the Higgs mechanism renders 
the non-Abelian gauge bosons massive. This remains true in the Higgs phase (broken 
phase) of the SM at high temperatures for small Higgs masses, as has been discussed 
carefully in the last years [Q-^. But due to the rising Higgs-plasma-mass at very high 
temperatures T > Tc ~ lOOGeV a new hot (unbroken) phase of the electroweak theory 
was argued to exist |0] in which according to perturbative calculations the gauge bosons 
are massless. A first order phase transition between the two phases could have important 
effects, in particular the possibility of a baryogenesis HI at this state has been discussed 
extensively j^. It now is practically excluded for the experimentally allowed range of 
rriH^QOGeV but variants of the standard model in particular the supersymmetric standard 
model are under discussion. 

If one wants to discuss the phase transition in detail ~ critical bubbles, sphaleron 
transitions, etc. - one needs the effective action (or at least the effective potential) of the 
theory. However, for vanishing or small Higgs vacuum expectation values the infrared 
problem comes up like in QCD and we expect non-perturbative effects. In the case of 
large m//>70GeV, where the phase transition starts to fade away, this non-perturbative 
regime might cover even the broken phase at the temperatures of the phase transition [10|. 

At high temperatures the non-zero Matsubara-modes are heavy and the IR-sensitive 
part of the theory is an effective three dimensional theory of the zero modes whose pa- 
rameters are obtained from the original theory by matching a set of amplitudes in 4 and 3 
dimensions [|^,^. The three dimensional action can be simplified further by integrating out 
the massive zero modes of the time component of the gauge field. In the high temperature 
expansion higher derivative terms are suppressed and neglecting electroweak mixing one 
ends up with the three dimensional SU(2)-Higgs model. An appropriate way to treat such 
an IR-sensitive theory where non-perturbative effects are expected is Monte Carlo simula- 
tions on the lattice like in QCD. This is simpler than in the latter case because the fermions 
have already been integrated out as non-zero Matsubara-modes. Such calculations have 
been performed by several groups |11-16|. They confirmed that there are confinement 
effects to be discussed in detail later on. An effective 3-dimensional gauge coupling ris- 
ing in the infrared was also argued for in the framework of exact renormalization group 
equations 0|. 

In this paper we will discuss the spatial correlation scales of the physical states, i.e. 
of the corresponding gauge invariant operators, in the 3-dimensional SU(2)-Higgs-model. 
The simple-minded perturbative picture would predict long range correlations because 
of the massless gauge bosons, but the neglect of IR-effects is obviously wrong. Based 
on 1-loop gap equations it was argued that one has just another Higgs-phase |18]. The 
predicted vector-boson-mass is much smaller than the one calculated in the gauge invariant 
lattice calculations (but curiously not much different from (Landau) gauge fixed lattice 
results [16|). In a recent letter we proposed a model for two-dimensional bound states 
of light constituents in the hot electroweak phase. The Green's functions of these bound 
states have been evaluated following a method which was developed by Simonov in 
4-dimensional QCD. In this paper we elaborate the model in detail. The evaluation of the 
Green's function is generalized and reorganized. The bound state correlation masses are 
first calculated analytically for a linear scalar-scalar-potential. Several modifications of 
the potential and their infiuences on the correlation masses are discussed. The parameters 
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of the potential are then fixed from lattice calculations of the Wegner- Wilson-loop. Finally 
we can compare our predictions for higher correlations with recent gauge invariant lattice 
calculations |13, T^. 



Chapter 2 introduces the effective 3-dimensional SU(2)-Higgs-model. Chapter 3 dis- 
cusses the possible bound states. The corresponding Green's functions are derived in 
chapter 4. They are evaluated in the next chapter. In chapter 6 we calculate the correla- 
tion masses for a linear potential, discuss modifications of this potential and their influence 
on the masses. The parameters of our model are fixed in chapter 7. Chapter 8 gives the 
comparison of our results with lattice data. The intercept of the potential is investigated 
in chapter 9. Finally we present our conclusions. 

2 The SU(2)-Higgs-model in three dimensions 

At high temperatures the SM can effectively be described by the three dimensional SU(2)- 
Higgs-model Q with the Lagrangian 

L = \ F^F^ + (A0)^(A</>) + ml 0t0 + A3(0V)' , A = a. - <734"y • (1) 

Even variants of the SM, e.g the minimal super symmetric standard model, can be described 



by this effective theory in a large part of the parameter space |21|. The squared mass 



m^, the quartic coupling A3 and the three dimensional gauge coupling depend on the 
parameters of the fundamental theory and on the temperature. The quartic coupling has 
the mass-dimension 1, while 53 has the mass-dimension 1/2. 

The non-perturbative aspects we are interested in are dominated by the gauge boson 
sector. In contrast to a four dimensional theory this sector has a natural mass scale 
given by the gauge coupling. It is therefore natural to express dimensioned quantities in 
units of powers of 53. The parameters of the model can then be represented by the two 
dimensionless quotients 



A3 = ^ and mi = ^ , (2) 



m^i — — 

93 93 



sometimes called x and y in the literature [y. The mass 771,3 depends in general on the 



renormalization scale. We work, however, with fixed values of 7773 and hence at a fixed 




scale. This is similar to lattice calculations where ml is fixed at the scale given by the 
lattice constant. 

In order to uncover the full symmetry of the SU(2)-Higgs-model we replace the Higgs- 
doublet (j){x) by a 2 X 2-matrix-field via 

= with * = ^0 12^2 + i^^T' , (3) 

where the r" are the Pauli matrices. Due to the identity Tr$^$ = 2 </)'''(/) the potential 
(eq. ([l|)) depends only on this trace. 

Expressing the gauge field by the usual matrix notation Ai = Af the action can 
be written as 

S3 = J d'x I ^Ti F.jFij + iTr(A*)^(A*) + Vi(^^TT^^^^ | . (4) 
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It is invariant under the transformation 

U{x) V , (5) 

U{x) Diix) ?7(x)t , (6) 

U{x) Fij{x) U{x)^ . (7) 

U{x) is an element of the local gauge group SU(2)g and depends in general on x, while 
y is a x-independent element of the isospin group SU(2)i. Hence the first index a of the 
matrix refers to the gauge group and the second index a to the global SU{2)j. 
has a first gauge index a and a second isospin index a. 



Di{x) 
Fij{x) 



3 Bound states of two scalars 

In the 3-dimensional SU(2)-Higgs-model like in all gauge theories in principle everything 
has to be expressed in terms of gauge invariant quantities. For the well known formulation 
of the Higgs mechanism for generating massive gauge bosons this just implies a rather 
trivial reformulation in terms of gauge invariant objects (unitary gauge). In the case of a 
confining theory such a gauge invariant description is standard and absolutely mandatory. 
Indeed some time ago the 4-dimensional SU(2)-Higgs model was a prominent model for 
describing an effective electroweak interaction of composite W, Higgs, quarks and leptons. 
But it is not clear just from a gauge invariant formulation whether one is in the Higgs 
or in the confinement phase. Contrary to the case of the original Abbot-Farhi-model for 
compositeness in electroweak interaction |22-24] we argue here that in the 3-dimensional 
case of the hot electroweak theory one really has a non-perturbative theory with genuine 
composite states in 2+1 dimensions. In this we differ from reference where a Higgs- 
mechanism also in this phase is advocated. 



Nonlocal operators The nonlocal operators corresponding to the bound state of two 
fundamental scalars $ are either 

Tr ^{xyT{x,x)^{x) (8) 

or 

Tt *(x)tT(x,x)*(x)r* . (9) 

The matrix T{x, x) with gauge-indices is the link-operator (also called gauge field trans- 
porter) defined by 

T{x,x) = rexp(^-ig3 £ Akdzk^ , (10) 

where V denotes the path-ordering operator along some given path between x and x. 
The operator of eq. (P) is an isospin-singlet, the operator of eq. is an isospin-triplet; 
both are gauge singlets. Since they are non-local all possible angular momentum states 
contribute. 

The local interpolating fields 

Tr *(x)"f*(x) and Tr *(x)^ *(x) (11) 

are included in these operators in the limit j; — > x. The index of the vector field corresponds 
to the spatial direction of the link operator. Note that Tr $(x)^$(x)r* vanishes while 
Tr ^{x)^Dk^{x) is a total derivaitve. 
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Selection rules The local limit shows that the isospin-singlet operator overlaps with 
scalar states while the isospin-triplet operator overlaps with vector states. The obvious 
question arises if there are scalar isospin-triplets or vector isospin-singlets as well. Let us 
first assume that the path between x and x is a straight line. In coordinate gauge the 
operator T{x,x) is then equal to the 11-matrix. The gauge-invariant operators and 
can be evaluated in this gauge. One gets 

Tr $(-r)t$(r)0/ with Oj e {12x2,t^} ■ (12) 

Here we have changed the coordinate system; the origin is now at the center of the bound 
state. 

The projection on parity-even and parity-odd states respectively is 

^Tr (*(-r)'''*(r) ± *(r)^ $(-r) ) O/ . (13) 
Expressing $(r) by real fields and as in eq. (^ the parity-even states become 

($o(-r)$o(r) + M-r)Mr)) Tr 12x2 O/ (14) 

while the parity-odd states are 

i($o(-r)$i(r) - $,(-r)$o(r-) - <^ki-r)Mr)ekli) Tr r' Oi . (15) 

One sees the parity-even isospin-triplets (with Oj = t^) and parity-odd isospin-singlets 
(with Oi = I2X2) vanish. Hence in particular no scalar isospin-triplets and no vector- 
isospin-singlets are allowed. 

If the path between x and x is not a straight line there may be an overlap with states 
which do no respect this selection rule. Nevertheless, there is no local interpolating field 
composed of two scalar fields alone which corresponds to these states in the limit x^x. 
The interpolating fields in that case would contain at least an additional gauge boson field 
Fij and thus correspond to hybrid states in QCD. The latter ones are several hundred MeV 
(order of magnitude of the "gluon constituent" mass) heavier than the corresponding pure 
quark-antiquark states. We assume that a similar mechanism holds also in our case. 

The derivation of the selection rules given above can be generalized to "blocked" op- 
erators which are used in lattice investigations. 



4 Green's functions and correlation masses 



In order to get some insight into the dynamical structure of the electroweak theory above 
the critical temperature and to compare with QCD in the confining phase we have pro- 
posed in reference |jl^ a bound state model with a potential whose parameters were de- 
termined by comparing with results of lattice simulations for Wegner- Wilson- loops. Since 
the fundamental scalars might be very light one has to treat the problem with relativistic 
kinematics. For that we adopted the method of Simonov |20| to our case. We reorganize 
and generalize our calculations of reference [|l9| in this paper. In the described model the 
correlation masses are determined from the exponential falloff of the Green's function. 
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The Green's function of a fundamental scalar We start from the Green's function 
G{x, y, A)aab/3 that transforms the matrix field ^{x)aa into ^{y)hj3- Neglecting the scalar 
self-interaction it has in the worldline formalism ||25tl the form 



G{x,y,A)aab(3 

dsexpf — mgs) / Vz 



(16) 



V exp 



dr ( \zI{t) + igsAkZk 



^ab 



a,l3 



where rrig is the squared mass of the fundamental scalar, s is the Schwinger proper time. 
Zk{T) is a path which connects x with y and is parameterized by r; r runs from to s. 
Zk^r) denotes the derivative of -Zfc(r) with respect to r. 

The path integral / Dz is performed over all paths from x to y. The measure of the 
integral is the exponential function of the wordline action 



^ exp ^- ^ dr Q (r) + igs A^Zk j 



(17) 



The first term in the exponent alone describes the free propagation of a scalar field. The 
second term includes the interaction with the gauge field A^ ; gs is the three dimensional 
gauge coupling. Scalar self-interactions are neglected (quenched approximation). The 
exponential function is ordered along the path; V denotes the path-ordering operator. The 
gauge field in eq. ( p^ ) is treated as a fixed background field. Later on we will quantize 
averaging over this field. (In principle there is one worldline action for every index pair 
(a,/3).) 

The path-ordered exponential function (eq. (p!7|)) is matrix valued. It is identical to 

^ ^ ' ^ Vexpf-igs J^drAkZk 



exp 



drzlir] 



a, 13 



exp(^-^^(iTifc(r)^ T{x,y)a,p , 



(18) 



where T{x,y) is the link-operator from eq. (10). 

Using this the Green's function (eq. ([IBD) can be written in Feynman-Schwinger-repre- 
sentation 

G{x,y,A)aabi3 = (isexp(-m|s) j Vz exp(^- ^ J^dr zI{t)J T{x,y)af3 Kb ■ 

(19) 

The path used for the evaluation of T{x, y) is the one integrated over in the path integral. 
All interactions of the scalar with the gauge field are expressed by this link-operator. 



The Green's function of a pair of scalars The Green's function G{x,x,y,y) con- 
necting the two isospin singlet operators Tr $(x)^ T{x, x) $(x) and Tr $(y)^ T{y, y) $(y) 
is built up by two fundamental Green's functions (eq. (p^)) and two link operators T 

G{x,x,y,y) (20) 
G{x,y,A)aab/3 T{y,y)i3^ G{y,x,A)b^aS T{x,x)sa 



-ig3 f Ak dzk 



2 J ds J ds exp (-m|(s + s)) JJvzVz 



10 Jo 
X exp 



V exp 



5 



The brackets {■ ■ ■) a symbolize the functional integration over the gauge field A. In 
principle this averaging contains all perturbative and non-perturbative effects. The last 
factor is the expectation value of the Wegner- Wilson-loop in a gauge field background. It 
results from four link-operators; two in the Green's functions of the fundamental scalars 
(eq. (p^) and two in the link-operators of the bound states (eq. |8|). The paths z{t) and 
z{f) determine the both sides of the Wegner- Wilson- loop. 
The Green's function of a triplet is correspondingly 

GT{x,x,y,y) (21) 

G{x,y,A)aabt3 T{y,y)i3^Tl^ G{y,X,A)cyd5 T{x,x)sa T^da)^ 

= S'^ G{x,x,y,y) (22) 

For a given i=j the Green's function of a triplet is identical to the Green's function of a 
singlet, but due to the selection rules discussed in section |3| the iso-singlet and iso-triplet 
states must have different angular momenta. 

The correlation mass The correlation mass of the bound states A4 describes the fall 
off of the Green's function at large distances Q 

G{x,x,y,y) oc exp(-7We) . (23) 

The extension of the states has to be small as compared to their distance 

|x — \y — y\ <C \x — y\, \x — y\ ^ Q . (24) 

An exact definition of M is given by 

M = - lim ^ln{G{x,x,y,y)) . (25) 

Prom the form of the Green's function G{x, x, y, y) it follows that M depends only on 
7713 and on the expectation value of the Wegner- Wilson-Loop. The other variables used 
in eq. ( pO| ) are integrated out. All possible angular momenta contribute to G{x,x,y,y). 
They have to be isolated later. 

5 The calculation of the correlation masses 



As in Simonovs original work [^0[ the correlation masses are calculated from eq. ( p5| ) 
by simplifying the Green's function of the bound state. The calculations are, however, 
reorganized and generalized. Furthermore we have in our (scalar) case no problems with 
chiral symmetry-breaking and spin-spin-interactions. In order to gain a manageable form 
we have to make two approximations. 

The parameterization of the paths In a first step the paths of both scalars are 
parameterized by a common parameter 7 via 

T = 7s f = 7s with < 7 < 1 . (26) 
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The kinetic terms of eq. ( pOD become 

The spatial vectors -2^(7) and ^^(7) are then replaced by a kind of "center of mass" 
and a relative coordinate 

ss / 1 1_ \ 

Rk = —z -Zk + -Zk Uk = Zk- Zk . (28) 

s+s \ s s J 

For Zk and Zk one gets 

s s 

Zk = Rk -\ — Uk Zk = Rk — Uk . (29) 

s+s s+s 



The path integrals / / DzDz in eq. (|20D transform into path integral over the new 
variables / / VRVu. and the kinetic terms (eq. (|27|)) become 



Mixed terms cancel. This parameterization is always possible and no approximation has 
been made. 



The classical path of the center of mass Now as in |2C| we make the crucial as- 
sumption that the classical path dominates the trajectory of the center of the bound state. 
This approximation is justified for heavy bound states. It is not obvious which scale has 
to be used as reference, therefore we can not decide here if the "large mass condition" is 



fulfilled; we will come back to this question in section 10. 

In this approximation the path integral / VR can be replaced by an appropriate pa- 
rameterization of Rk{'y)- We choose the xs-axes as the direction of the path 

R^ = R2 = R^ = = 'd ^ (31) 



where is the distance between the two bound states introduced in eq. (p4[). The param- 
eter I? will be interpreted in connection with eq. (|37|). 
The Green's function of the bound state becomes 

G{x,x,y,y) (32) 
oc J ds J ds exp^— m|(s + s)j 
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The modified area law The Wegner- Wilson-Loop is unchanged by the former trans- 
formations. Motivated by lattice results we propose the ansatz 



V exp -ig3 



Ak dzk 

x,x,y,y 



OC exp 



(- e J^'dj V (v'nf(7)+ni(7)) ) (33) 



for the gauge field averaged expression. This replacement generalizes the evaluation of 
the Wegner- Wilson- loop in the references [20| and [l^. As we will see in section 3.2, 
effects important at small distances r can be taken into account (cf. sect. ^^) with this 
generalization. 

Note that the averaging over the gauge field is included in the modified area law 
(eq. (33)). The non-perturbative effects originate in the gauge boson sector and are repre- 
sented by this averaging. We thus assume that all non-perturbative effects on the bound 
state masses can be expressed by a scalar-scalar-potential V{r). 

Eq. (|3^) evaluated for rectangular Wegner- Wilson-loops with infinite length T leads 
to the two-dimensional potential 

V{r) = - lim ^ln{W{r,T)) . (34) 



T 



Here W{r, T) denotes the expectation value of the rectangular loop with length T and 
width r. In section ^ we use this definition to determine V{r) from lattice data. 

Using the replacement of eq. (p3[) the ii3-path integral in the Green's function of the 



bound state (eq. (|32|)) becomes independent of 0. Therefore it does not contribute to the 
correlation mass. Since the Green's function is only needed up to a constant we disregard 
this contribution. The two remaining coordinates are combined to a vector u = {ui,U2) ■ 



Evaluation of the Green's function In the next step we find a Hamiltonian for a 
corresponding Schrodinger type equation whose Green's function is the one of eq. ( ^2|) 
with the replacement of eq. (^). The latter should bahave as G oc exp(— Mj0) in the 
channel i. This has two advantages: First the separation of states with different angular 
momentum is very simple and second the solution of the resulting differential equations is 
much simpler than performing the path integrals. 

For this purpose it is advisable to introduce new variables. First 7 is replaced by the 
parameter '& from eq. (]3l|). Then the Schwinger proper time variables are substituted by 

e _ e . e 

fj, = — fJ- = — A^ = = 777 =7 • ("j5) 



2s 

The Green's function (eq. 

G{x,x,y,y) 



2s 

|)) becomes 



+ 



2(s+s) 



(36) 



oc 



e 
2" 



I 1 

- + - 

II II 



+ (m + m) 



Du exp{—B) 



with 



B 



e 



/ du{'&) V 
2 \ dd~) 



+ vm 



(37) 



If we look first to the -u-path integral the formal analogy to a two dimensional system 
is obvious: a particle with mass jl at position u and time "i? moves in the potential V {\u\)- 
The Euclidean action is B; the time interval runs from to 0. Hence, the component 
?? = i?3 of the center of mass coordinate takes over the role of the quasi time. 
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The Schrodinger equation of the two dimensional problem. The time indepen- 
dent Schrodinger equation corresponding to the two dimensional problem is 

H'^iu) = e*(n) (38) 

where e is the energy eigenvalue and the Hamilton operator is given by 

Due to the rotational symmetry of H we can separate the radial from the angular variables 
and the problem simplifies to the solution of the radial equation. Using r = |u| one gets 

- ^[7^ + 2]'f'nlir) + V{r) Mr) = e^i Mr) , (40) 

2/i \ ar^ r or j 

where ?7- = 1, 2, . . . is the radial and / = 0, ±1, ±2, ... is the angular momentum quantum 
number. 

If V{r) is less singular than 1/r at r = 0, the boundary conditions are 
lim Tpnlir) = and ipnlir) = at r = . (41) 

r — *oo 

For a given potential V{r) and fixed quantum numbers n and / it is possible to solve the 
radial equation and to calculate the eigenvalue Cni- In section ^ we derive the results for 
some potentials of interest. 



The Green's function of a particular state The solution of the Schrodinger equation 
(^8|) allows us to express the n-dependent part of the Green's function by 

Jvuexp{-B) = rnMO)) Mm) eM-<^niQ) (42) 

n,l 

Thus the contribution of a particular state with quantum numbers n and / to the 0- 
dependence of the Green's function (eq. (|3^) ) is given by 



Gni{x,x,y,y) (43) 

^ lo lo (2^) ^^^(~^ "^Km"^^) + 2e„i(2/i) ^ . 

The saddle point method The fi- and the /i-integrals are evaluated with the saddle 
point method. For large values of Q both integrals are dominated by the exponential 
function. In the limiting case of 0^oo only the maximum of the exponent with respect 
to fi and /Z contributes. 

The calculation may by abbreviated using the fact that the exponent is symmetric in 
fi and p.. This is due to the fact that eq. ( pOD is symmetric in s and s. This symmetry in 
turn originates from the identity of the masses of both fundamental scalars. The saddle 
point equations for fj, and fi are hence equivalent. The minimum of both is unique and 
therefore the same. It is possible to identify both parameters. Using /i = /i = 2/i (cf. eq. 
^ ) we arrive at the Green's function 

Gni{x,x,y,y) oc exp{- Mniil^) &) (44) 
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with the correlation mass 



where fj, fulfills the saddle point equation 







(45) 



(46) 



Thus the calculation of the correlation mass from the Green's function of the bound 
state (eq. (pO|) ) is reduced to the solution of a differential equation (eq. (^0|)) and a mini- 



mization problem (eq. (46)). 

6 The correlation masses 
6.1 Linear potential 

The solutions of the Schrodinger equation (eq. (p8|)) and the corresponding eigenvalues 
Cni depend on the scalar-scalar-potential. In this subsection we investigate the simple and 
important case of a linear potential 

V{r) = ar . (47) 
The radial equation Substituting r in eq. (|40| ) by the dimensionless variable 



(^CJ) 



1/3, 



results in the eigenvalue equation 



P dp 



+ 



+ p]lpnl{p) = anllpnlip) 



(48) 



(49) 



The boundary conditions (eq. (|4l|)) are unchanged. 

The numerical solution of this differential equation does not cause any difficulties. The 
eigenfunction with the lowest quantum numbers are displayed in figure ||. The correspond- 
ing eigenvalues ani are: 





1=0 1=1 1=2 


n=l 
n=2 


1.74 2.87 3.82 
3.67 4.49 5.26 



The eigenvalues of interest €ni{p) can by calculated from the dimensionless eigenvalues 
ani of eq. (H) via 

^2/3 



P 



1/3 



a-ni 



(50) 



Therefore the /x-dependence of €ni{p) is a power law in the case of a linear potential. 
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Figure 1: The solutions of the radial equation (|49|) for the quantum numbers n = 1, 2 



and / = 0, 1, 2. The Green's function of the bound states (36) can be expressed be these 
solutions. The notation s-, p- and d-wave corresponds to different angular momentum 
quantum numbers /, as in atomic physics. 

The saddle point equation Using e„;(^) from eq. (|50| ) the saddle point equation ( ^6|) 
becomes 

- 4 + 1 - ^ M-'/V2/3a„, = . (51) 

It is solved by 

2 

H = zianuhf'"^ Vo with b = ^ (52) 

a 

where z(a, b) is a positive and real solution of the cubic equation 

- ^az - b = . (53) 
The solution corresponding to the minimum of Ai (fi) is 

2I/3 1/3 1 

z{a,b) = —a (27b + V72962 - AaA " + — — (276 + ^729^2 - 4aA . (54) 

It is unique if 6>0; for —-^a'^^'^<b<0 it corresponds to a local minimum of 7W. If 6 is 
even smaller there is no positive and real solution anymore. We therefore treat only the 
case that ml and hence b is positive. 

The correlation mass (eq. (|45|)) evaluates to 



Mni ( m'i\^^'^ ( ■mV\ '^^'^ ml , , 

= 4z a„z,^ - 2z a„i,-3 -3 , (55) 



with the function z given in eq. (^). It depends only on the dimensionless eigenvalue 



anil on the string tension a and on the squared mass of the fundamental scalar. The 
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2.0 ' • ' • ' • ' • ' • ' 

0.0 0.2 0.4 0.6 0.8 1.0 

ml /a 

Figure 2: The correlation mass M.ni of the bound states versus the squared mass of the 
fundamental scalar ml for a linear potential. We find a spectrum of higher excitations. 
The string tension a sets the scale of the bound state masses. In the case of a linear scalar- 
scalar-potential the dependence of the bound state mass on the mass of the fundamental 
scalar can be given analytically. 

result scales with a. In figure || we plot the correlation masses of the bound states with 
the lowest quantum numbers versus the squared mass of the fundamental scalar. 

Our results show a hierarchy of bound states. The lowest state corresponds to the 
Is- wave (n = 1, / = 0). Due to the selection rules derived in section ^ it has to be an 
isospin-singlet. The Ip-state {n = l, 1 = 1) is heavier, it is a vector-isospin-triplet. A dense 
spectrum of higher excitations follows. Note that the string tension determines the typical 
scale of the bound state masses. It is a genuine non-perturbative quantity. 

For the case that ml<a our treatment is by no means an analogue of the non- 
relativistic quark model . This can be seen by comparing eq. (|40| ) with eq. (^). From 
( p5[) we would deduce a constituent mass of the elementary scalars 

fi ml 

"Iconstit = + ^^^) 

whereas from the kinetic terms we would deduce ?7iconstit = fJ- ■ For ml>a there is no 
contradiction since /i goes to ml in that case. Furthermore the constituent mass shows a 
marked dependence on the state for ml<a . 

The mass parameter /x varies from 0.2Agl to OAgl going from Is to 2p corresponding 
to a quasi static constituent mass of half that value at ml = 



^Thus it is about half of the constituent mass introduced in a very recent paper by Buchmiiller and 
Philipsen [Q. Note that we also have a binding energy e„;. The constituent W-boson in does not enter 
directly the mass formula in our dynamical approach. However it corresponds to the inverse correlation 
length of the gauge field strength and hence determines the binding energy via the string tension (eq. (^) 
and appendix 
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a) b) c) 




Figure 3: Different modifications of tlie linear potential, a) the influence of screening, b) 
the lowering of the potential due to non-perturbative effects, c) a perturbative contribution. 



6.2 A nonlinear potential 

The procedure presented above can be extended in a straight forward way to a general 



potential |26|. For the problem under consideration several modifications of the linear 
potential are expected to occur. 

As in QCD we do expect screening of the potential by spontaneous creation of a pair 
of two fundamental scalars. The potential would look like the full line in figure |3|a. The 
dashed line shows the corresponding linear potential. So far lattice results do not indicate 
any screening and the low lying bound states should not be affected by the screening 
anyway. We therefore neglect this effect. 

From lattice calculations the potential can only be determined up to an additional 
constant. In the stochastic vacuum model, outlined in appendix |^, the linear rise of the 
potential sets in only at distances r which are larger than the correlation length of two 
gauge fields in the vacuum. The shape of the potential is given in figure ^3 (full line). 
If the wave functions obtained from the Schrodinger equation are not dominated by the 
region of small r, the potential can be approximated by a linear potential with a negative 
intercept 

V{r) = Vo + ar with Vb < (57) 

(dotted line). The size of the intercept will be discussed in section ^. 

It is easy to see that the correlation masses corresponding to this potential are 

Mni = M'^ + Vo , (58) 

where are the masses calculated from the linear potential. Due to the lowering of 

the potential the correlation masses of the bound states are lowered as well. 

The difference between the real potential and the approximation (eq. (p7[)) may, if it 
is small, be taken into account as a perturbation (see below). 

At small distances r the potential should be dominated by the perturbative gauge 
boson exchange. The lowest order is calculated in appendix This leads to the two 
dimensional Coulomb potential 

= ^glHA\u\) . (59) 
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The constant A cannot be fixed due to the logarithmic divergence. 

If the exchanged particle acquires the mass m we obtain the two dimensional Yukawa 
potential 

^y(I^I) = -^9lKo{m\u\) , (60) 

where Kq is the modified Bessel function. Note that both potentials merge at small 
distances \u\ = r<^m~^ with a suitable choice of A. At these distances the mass does not 
make any difference. 

If the perturbative contribution to the potential vanishes in the limit r— >oo as VY(r) 
does, the linear potential is not changed here. This situation is sketched in figure |^. It 
turns out that in this case the corrections to the correlation masses Aini can be calculated 
perturbatively. One gets 

Mnl = Mlf + dMnl (61) 

with 

roc 

6Mni = 27r / drrVY{r)iP^i{rf . (62) 

JO 

Since the Yukawa potential is negative, the masses are lowered. 



7 The potential on the lattice 

In order to compare the correlation masses of the bound states with lattice data, as will 
be done in the next section, we have to determine V{r). In this section we analyze lattice 
data of Wegner- Wilson-loops and fix the parameters of the potential using a suitable fit. 

The symmetric electroweak phase of the three dimensional SU(2)-Higgs-model has 
been investigated by several groups on the lattice |Tl|-^. Nevertheless, only the authors 
of reference |ll|-||] calculated the We gner- Wilson-loops. The data have been reanalyzed 



by us 1 26 1 and will be given in a form which is appropriate for the purpose of calculating 
the bound state masses. 

The lattice data are expressed in terms of the parameters /Jg, and (in. For j3G 
there is the simple relation 

Pg = , (63) 

where ai is the lattice constant. The data we use are evaluated at /?g = 12 and hence 
correspond to = ^ g^'^ . The more complicated relation between the other lattice 
parameters and the (normalized) continuum parameters m\ = m'^/g^ and A3 = A3/53 
(eq. (§)) are given in reference |2^. The data of correspond to A3 = 0.0239 and 
mi = 0.73; the A3 = 0.0957 results cover the range m| = -0.022 to 0.524. 



The Wegner- Wilson-loops on the lattice are rectangular and have extensions from 2 
to A^ = 15resp. 24 lattice units. The extrapolation to infinite length is very secure, since 
the relevant quantity — In VF(r, T) (cf. eq. jsj) rises practically linear in T for T> 8a. We 
found that the lattice data can be reproduced by the three parameter fit of the form 

V{r) = C - —glKoimr) + ar (64) 
8 vr 

The fitted constant C is without any physical significance since it depends on the lattice 
renormalization procedure. An effective intercept C = Vq like in eq. ( [57| ) due to non- 
perturbative effects will be discussed in section ^. In table || we give the results for the 
relevant parameters m and a for the available values of A3 and m'^. 
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A3 


9 

mj 


N 




/ A 


0.0957 


0.5254 


15 


1.045 


0.1370 


0.0957 


0.3699 


15 


1.105 


0.1390 


0.0957 


0.2153 


15 


1.029 


0.1345 


0.0957 


0.1229 


15 


1.083 


0.1364 


0.0957 


0.0615 


15 


1.046 


0.1335 


0.0957 


0.0309 


15 


0.999 


0.1296 


0.0957 


0.0003 


15 


0.886 


0.1208 


0.0957 


-0.0058 


24 


0.770 


0.1135 


0.0957 


-0.0220 


24 


0.411 


0.0876 


0.0957 


0.073 




1.058 


0.1345 


0.0239 


0.073 


15 


1.051 


0.1326 



Table 1: The fitted parameter of the function (eq. |6J). The values at A3 = 0.0957 and 
ml = 0.073 are interpolated. is the number of lattice points, m is the effective gauge 
boson mass and a is the string tension. 

The contribution -^gk^KQ{mr) indicates that the coloured objects exchanged between 
the scalars have an effective mass of about 1^1 . It is suggestive to identify it with some 
magnetic mass of the gauge boson. The mass m of our fit is, however, larger than the 
magnetic mass obtained in lattice simulations in Landau gauge |16] and also larger than 
the one predicted in reference ||l^ from gap equations. 

Comparing the potential parameters at A3 = 0.0239 with the ones at A3 = 0.0957 
(interpolated) one sees practically no effect of the quartic coupling. This can be understood 
from the fact that the exchange of a scalar is suppressed by a factor A3 in comparison to 
the gauge boson exchange. 

The shape of the scalar-scalar-potential depends at > 0.0615 only marginally on 
Mg. In this mass range V[r^ is nearly exclusively determined by the gauge boson sector. 
This behavior is expected for large scalar masses, since the scalars decouple here. From 
the data of table |l] we conclude that the mass-limit from which on the influence of the 
scalars on l/(r) is negligible is small compared to the other mass scales of the problem. 
Above this limit the string tension a should be identical to the one of the pure SU(2)- 



Yang-Mills-Theory. Indeed, Teper |28| finds for the latter one roug hly a 0.137^1 at 
= 12 based on the correlation of Polyakov- loops. This value is similar to the string 
tensions of tabel He extrapolates his results to /?g = c>o in order to gain a continuum 
value of the strig tension. We stay, however, with the potential at /^g = 12 , since we 
compare with lattice data at the same value in the next chapter. 

The data point ml = —0.022 is at the critical temperature. The fact that the critical 



''3 

t2 \ 



ml is negative is in agreement with other lattice investigations [14|. 
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As 


ml 


*^ is 


^ ^ ZS 


^ *■ I J^p 




g'i 


0.0957 


0.5254 


1.93 


2.48 


2.27 


2.71 


2.53 


0.0957 


3699 


1.73 


2.31 


2.09 


2.55 


2.36 


0.0957 


0.2153 


1.48 


2.08 


1.85 


2.33 


2.14 


0.0957 


0.1229 


1.31 


1.95 


1.70 


2.21 


2.00 


0.0957 


0.0615 


1.15 


1.82 


1.56 


2.09 


1.88 


0.0957 


0.0309 


1.05 


1.74 


1.47 


2.01 


1.79 


0.0957 


0.0003 


0.91 


1.61 


1.34 


1.88 


1.67 


0.0239 


0.073 


1.18 


1.84 


1.58 


2.10 


1.89 


0.0239 


0.089 


1.22 


1.87 


1.62 


2.13 


1.92 



Table 2: The masses of the bound states up to an arbitrary additional constant C . They 



have been calculated using the potential eq. (|64D with the parameters of table Q and C = 



8 Comparison of the correlation masses with lattice calcu- 
lations 

The correlation masses in our model We calculated the correlation masses Mni of 
the bound states with the method presented above. We based this calculations on the 
potential of eq. (|6^ ) with m and a from table |l| and the arbitrary choise C = . The effect 
of the term —-^glKQ{rnr) is very small. It lowers the masses by less than 2% and can 
safely be treated perturbatively. This is due to the large extensions of the bound states 
as compared to 1/m. Our results for the correlation masses are listed in table ||. 

In lattice calculations the correlations of pairs of the operators (eq. (^)) respectively 
(eq. d^) are investigated. The relative orientation can be varied. Using special combina- 
tions of these pairs it is possible to extract single angular momentum states. The method 
of blocking operators and the diagonalization of mass matrices are used to suppress the 
mixing of different states. In this way one gets the correlations masses corresponding to 
definite spin and isospin quantum numbers. 

The quantum numbers of the Is-state are the same as those of the Higgs-particle; the 
quantum numbers of the Ip-state are identical with those of the T^-boson. Therefore the 
operators used to investigate the Is- and the Ip-state are the same operators which are 
used in the broken phase to investigate the Higgs- and VF-boson, respectively. Therefore 
the Is-mass is sometimes called Higgs-mass, the Ip-mass is called Pl^-mass in the literature. 

Philipsen, Teper and Wittig |l5| calculated the correlation masses of the bound states 



at As = 0.0239 and ml = 0.089 on the lattice. They find a spectrum of higher excitations. 
The masses of the low lying states are determined with high accuracy 

TWPTW ^ (0.839 ±0.015)51 , 
A^r^w ^ (1.27 ± 0.06) 5i , 
Mir = (1.47 ±0.04)^1 . 

Since the Wegner- Wilson-loops have not been measured by Philipsen et al., we have to 
use the potentials analyzed above to compare with our results. We do not have any data 
at the values of A3 and ml used in reference [15|. In view of the small 7713-dependence of 
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Figure 4: The correlation mass of the Is-state vs. the squared mass of the fundamental 
scalar. The points with error-bars are the lattice results of reference [jl^.The full line 
connects the results of our model. Both data sets are normalized at mi = . 



the potential (cf. tab. ||) we use instead the potential corresponding to A3 = 0.0239 and 
ml = 0.073 to calculate the correlation masses. The results are given in the last line of 
table |. 

Since unfortunately the constant C of the potential cannot be deduced from the lattice 
data, we cannot compare the masses directly, but only the mass differences, as given in 
the table 





{Mip-Mis)g3^ 




our model 
Philipsen et al. 


0.40 
0.43 ± 0.06 


0.65 
0.63 ±0.05 



Within the errors the data predicted by our model agree very well with those from the 
lattice. 



Giirtler et al. [13| have calculated the Is- masses as well as an upper bound on the 
Ip-masses at A3 = 0.0957. There results for the Is-masses for different scalar masses ml 
are plotted in figure Q in comparison with our results. Both data sets are normalized 
at ml = . One sees that the dependence of the bound state mass on the mass of the 
fundamental scalar is well described by our model. The main contribution to the variation 
of Ails is due to the explicit ml-dependence of the Green's function (eq. (|20|)), while the 



modification of the potential with m^ is not significant. 



9 The intercept 



The constant C in eq. (^4[) can not be fixed from lattice data of the Wegner- Wilson- 
Loop. One could choose it to give good agreement with the results of reference [13|, 



resulting in C = —OAAg^ . Similarly one could choose the constant to give the best 
possible agreement with the masses from reference [fH]. This would lead to the slightly 
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larger value C = —0.38 . The difference between these two values can in our opinion not 
be attributed to the different values of the quartic coupling. As we showed in section |^ the 
latter has only a small influence on the potential. It is rather due to a small discrepancy 
between the two lattice investigations, which is, however, within the statistical error. 

In order to get some theoretical understanding of the origin of the intercept it seems to 
be necessary to investigate the vacuum structure of the theory. In the bound state model 
we have considered the binding forces between the fundamental scalars. The parameters 
were taken from lattice calculations of the Wegner- Wilson loops. The resulting spectra 
strongly support a picture of the effective 3-dimensional electroweak theory above the 
critical temperature which is very similar to QCD with light scalar quarks, the linear 
confinement playing an essential role. Though there is no analytic proof for confinement 
in (4-dimensional) QCD there exists a simple model which yields confinement for non- 
Abelian gauge theories in a very natural way, the model of the stochastic vacuum 

mil- 



In appendix |^ we give the essential features of this model for a three-dimensional theory. 
The model yields an asymptotically linear potential and relates the vacuum expectation 
value of the gauge boson fields ( g^FF ) and their correlation lengths a to the string tension. 
For small inter-quark separations the potential is quadratic and becomes linear only at a 
distance of a few correlation lengths a and thus corresponds to the situation depicted in 
figure 3b. If we denote by D{z'^) the (rotationally invariant) correlation function of the 
gluon fields in the vacuum (see eq. (|79|)) the string tension is given by: 

whereas the effective intercept C = Vq (see figure 3b) is given by: 



oo 



Vo = -illZll dzz'D{z') (66) 



6Nc Jo 

From these equations we can deduce: 

Vo = -Kaa (67) 



where a is defined by 



/■oo 

/ dzD{z^) (68) 
Jo 



and K is of the order 1, the numerical value depends on the form of D{z'^). Unfortunately 
neither the condensate ( FF ) nor the correlation length have been calculated on the 
lattice so far, so we apply here a simplifying argument in order to get an idea of the order 
of magnitude of Vq. 

If a is the correlation length between two field strengths F we expect the correlation 
between the product of two field strengths F^ to be of order a/2.0 Since the product 
interpolates a glue-ball (VF-ball) we expect 2/ a to be near the glue-ball (VF-ball) mass. 
In QCD this seems to be fulfilled qualitatively, with 2/a~1.5 GeV. The I^-ball mass has 
been obtained in lattice calculations [|l5| and found to be 771(5 = 1.60 ± 0.04^3 . Assuming 
the correlation length to be a = 2/mG we obtain together with the string tension of the 
lattice calculations (cf. tab. |l|) the value Vq = —0.17K g"^ which is indeed of the same order 

■^Taking o/4 (corresponding to 4 gluons in the glueball |Q) would improve the agreement discussed 
below. We thank O. Philipsen for a discussion of this point. 
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of magnitude as the constant C which we need to fix the absolute values of the bound-state 
masses. It would be extremely interesting to determine the correlator ( F(x)T(x, 0)F{0) ) 
on the lattice in order to see if considerations similar as in QCD lead to confinement in 
this 3-dimensional theory. 



10 Discussion and conclusions 

The comparison with the lattice data in the last section shows that the correlation masses 
of physical (gauge invariant) objects in the hot electroweak phase are described well by the 
proposed model. It turns out that it is similar to QCD but without the problem of chiral 
symmetry breaking and spin-spin- interactions. The fact that the dense spectrum of higher 



excitations predicted by us in reference [19| has been confirmed by lattice calculation is 



in our opinion a success of the model on its own. The numerical agreement of the mass- 
differences is much better than it could be expected in view of the approximations made. 
Even the m^-dependence of the Is-mass is explained by the model. In consideration of 
this fact we will discuss the assumptions of our calculations again. 

The Green's function of the fundamental scalar (eq. (|T^)) does not take into account 
the interactions with other scalars (quenched approximation). Nevertheless, the influence 
of the dynamical scalars is taken into account. The fact that the scalar-scalar-potential 
measured on the lattice is not the same for all data sets is due the dynamical scalars. 
The As-dependence is small, as far as we can conclude this from the existing data. The 
modification of the potential as function of the mass is only large for very small values 
of this mass. The direct perturbative exchange of a scalar can be taken into account 
like the perturbative gauge boson contribution. It is, however, suppressed by a factor A3 
compared to the latter one. The scalar contributions to the scalar self-energy is neglected 
by the quenched approximation as well. All these terms are suppressed by at least a 
factor A3 when compared to those contributions taken into account. This approximation 
is therefore well-founded. Only the influence of the dynamical scalars of the scalar-scalar- 
potential for very small squared masses is of numerical relevance; this effect is included 
in our calculations. 

In the calculation of the correlation masses from the Green's function (eq. (pO|)) in 
section || we made two approximations: the choice of the classical path for the center of 
mass trajectory and the modified area law. 

By the restriction to a straight and steady movement of the center of mass R{j) 
fiuctuations of the bound state as a unit are excluded. This approximation is well-founded 
for heavy bound states. As we have seen, the typical scale of the problem is the string 
tension. The masses of the bound states are A4 > 3^/a (cf. fig. ^) and hence larger than 
this scale. It might, nevertheless, be interesting to check this approximation or to calculate 
the next order corrections. Note that the relative movement of both fundamental scalars 
is treated fully relativistically in spite of the restriction to the center of mass path. 

In this context the obvious question arises why it is possible to reduce the treatment 
of the relative movement to a Schrodinger equation. Here it is important to note that it 
is only an analogous non-relativistic problem which is described by this equation. The 
parameters ^ and ft then have the meaning of constituent masses. They are averaged 
by the Schwinger-proper-time integral later on. The evaluation of this integral with the 
saddle point method is exact in the limit Q ^ 00 investigated by us. Neither the non- 
relativistic treatment of the analog problem nor the saddle point method is based on any 
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approximations. One may, at most, ask if the limit Q to infinity can be extrapolated by 
the lattice investigations we are comparing with. 

The modified area law is included in our calculations of the correlation masses. The 
correlation with the usual area law was already discussed in section ^. An extension to 
a potential which depends on all three components of u does not cause any problem, is, 
however, not adequate since the connection with the scalar-scalar-potential calculated on 
the lattice would get lost. 

Finally we neglected the effect of the screening. The Wegner- Wilson-loops calculated 
on the lattice show no screening of the potential up to distances of about 8g^^. In figure 
^ we plotted V{r) in comparison with the bound state wave functions. One sees that the 
wave functions corresponding to the lowest quantum numbers have extensions between 
8g^'^ and lOg^'^ . To exclude screening effects at all the linear part of the potential has to 
be confirmed at distances which are slightly larger than those of todays data. It is, on the 
other hand, not expected that the lowest excitations are influenced by screening effects. 

The very good agreement of the correlation masses calculated by us with the lattice 
data justifies our approximations a posteriori. As long as the errors of the lattice data do 
not get smaller than today we see no need to go along without these assumptions. 

A more interesting question is the calculation of the scalar-scalar-potential. The Bessel- 
function contribution (eq. (|6^)) seems to originate in the exchange of a colored massive 
particle or quasi-particle. The mass is of order g^; the nature of the exchanged particle is 
not yet revealed. As explained in appendix ^ the string tension and the intercept can in 
principle be calculated with the model of the stochastic vacuum. The gluon correlator is 
an important input function for this model, which needs confirmation from the lattice. 

Another interesting point, which is not fully clarified yet, is the treatment of negative 
values of 777,3. It is not attractive to have a fundamental scalar with negative squared mass. 
One could try to use the perturbative mass calculated from the second derivative of the 
effective potential at (/? = instead. The latter one diverges, however, due to the ip"^ ln((/?)- 
term of the potential. This infrared-divergence shows the breakdown of the perturbation 
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theory in the symmetric phase. A much more speculative possibihty is the generation of 
a scalar mass by a gluon-condensate, similar to the generation of a gauge boson mass by 
a scalar condensate in the broken phase. A gluon condensate due to the instability of the 
naive vacuum for small Higgs vev's should also be very important for the generation of 
the non-perturbative part of the electroweak effective potential. 

An important outcome of this investigation is that the hot electroweak phase is strongly 
interacting like a pure 3-dimensional Yang-Mills-theory. The spectrum of gauge invariant 
states is, however, totally different from the latter one. The existence of low-lying Is- and 
Ip-states is only due to the fundamental scalars; the gauge invariant spin 0- and spin 1- 
operators of eq. ( pTD do not exist in pure Yang-Mills-Theories. The gauge bosons in the 
pure gauge theory are only defined in a definite gauge while the gauge-invariant spin 1- 
object is a W-ball. Indeed there is the interesting observation from lattice calculations 
[p!^ , ^ that glueballs and Higgs states are practically decoupled. The string tension is 
universal in a given confining theory and thus the former remark fits very well together 
with the fact mentioned above that the string tension responsible for ^-binding is nearly 
the one of pure Yang -Mills theory in lattice calculations. 

Acknowledgement 

We would like to thank M. Ilgenfritz, M. Laine, O. Philipsen and M. Rueter for useful 
discussions. 



A Perturbative contributions to the potential 

Using the cumulant expansion |32| the expectation value of the rectangular Wegner- 
Wilson-loop can be expressed as 

1 „„.,.,.,».„) . 

(69) 

To lowest order in this expansion one has 

lnTy(r,r) = jdx,jdy,{{A-{x)A-{y))) . (70) 

The perturbative contribution is now calculated from the exchange of a perturbative gauge 
boson between the two long sides of the loop. With x and y as two-dimensional vectors 
and xi and yi the corresponding components running along the long sides one finds 

T/2 -T/2 

hm i In W{r,T) = - ^ lim i / dx^ f dyi {{A'i{x^,x)AUyi,y))) (71) 

-T/2 T/2 

„2 j'oo 

= f / dy^{{AUO,x)AUyi,y))) , (72) 

where we have used translational invariance in the last line. 

The two-point-cumulant is now replaced by the perturbative propagator in Feynman- 
gauge via 

J2^ WTT^ '^'^"'^^ ■ ^^'^ 
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We split the momentum coordinates in a first component ki and a two dimensional vector 
k. The perturbative part of the potential is 

VpeAy-^) = i rdyif4^^-^. e^i'^^y^^'^^y-^^^) (74) 



CO 



(27r)3 kf + k'^ + m^ 



(27r)2 ^2 + ^2 

It is thus the two-dimensional Fourier-transformation of the perturbative propagator in 
momentum space with fei = 0. The derivation holds for m = as well. 



B The model of the stochastic vacuum in three dimensions 

We introduced a model of bound states in the hot electroweak phase and showed how to 
calculate the correlation masses within this model. The knowledge of the scalar-scalar- 
potential V{r) is an important requirement of this calculation. The origin of the asymp- 
totically linear potential itself has not been explained so far. This can be attempted by 
adopting the model of the stochastic vacuum, which was originally proposed by Dosch and 
Simonov [p9| , |30| for QCD, to the three dimensional SU(2)-Higgs-model. 

The basic ingredient of the model of the stochastic vacuum is the assumption that the 
complicated contributions of non-perturbative field configurations can be approximated 
by a simple stochastic process. In that way already the assumption that this process has a 
convergent linked cluster expansion leads very naturally to linear confinement. Making the 
more restrictive approximation of a Gaussian stochastic process leads to a very predictive 
model since now in principle the full non-perturbative contribution is approximated by 
a single correlator. This correlator can be used to determine observables as the static 
quark-antiquark potential in QCD, soft high energy cross sections and others. It can also 
be compared with lattice calculations and all results turned out to be very satisfactory. 

In the following we shortly exhibit basic features of the model adopted to three di- 
mensions and SU(A'c) as gauge group with Nc = 2 , the case considered here. In three 
dimensions we can introduce besides the field tensor Ffi the vector 

Fk = ^r,kF^^ (76) 
which in an Abelian theory satisfies the Bianchi identity 

dkn = Q . (77) 

In order to form gauge invariant correlators of the field strength tensors we have to trans- 
port the color content of all fields at point x to a fixed reference point y. This is done by 
the gauge field transporter T{x,y) (eq. ([lO|)) along a straight line from x to y. 



i^5(x,y)y = T{x,yr'F,-j{x)^T{x,y) . (78) 

The correlator {glF^-{x,y)F^j^{x' ,y))A depends in general on the reference point y. 
We make the assumption that it can be approximated by an expression depending only 
on the difference z = x — x' . 



{F-{x,y)F',i{x',y))A = g(jv|-l) ^^^^ 



K {SikSji - 6u6jk) D{z^) (79) 
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Here {FF) is the gluon condensate 



(PF) = EE(^^-(0)^"(0)U (80) 

ij a 

and D and Di are correlation functions which are supposed to fall off with a characteristic 
correlation length a. 

In an Abelian gauge (without monopoles) the Bianchi identity (eq. (|77|)) forces the 
first contribution to vanish, i.e. k = , but in a non- Abelian theory we have no longer the 
Bianchi identity and hence there is no reason that the correlator D{z'^) vanishes. Lattice 



results [33| showed k « 0.74 in QCD. No lattice results are available for the 3-dimensional 
correlator. 

From the correlator one can obtain the area law by first transforming the line integral 
over the potential Af into a surface integral over the field strength F^j by means of the 
non-Abelian Stokes theorem and then applying the cluster expansion. For subtleties due 
to the path ordering we refer to reference . 

For a rectangular Wegner- Wilson-loop of extension T in zg-direction and r in zi- 
direction we obtain: 



In H^(r,T) = IniV, 



c 



g'jFF) 



i-T i-r 

/ dzo / dzi K {rT — rzo — Tzi + zqzi) D{z'^) 
Jo Jo 



+ ( ^Tzi + ^rzo - ziZQ ) Di[z^) 



(81) 



For small values of r and T the non-constant term in In W behaves as —-^^j^g^^FF) r^T^, 

independent of the value of k, whereas for large values of r and T only the correlator D{z'^) 
gives an area law. 



The static potential is obtained from W{r,T) by the limit of eq. (34). It increases 
quadratically in r for r<^a and linearly for r^a (cf. fig. ^b). Asymptotically, i.e. for 
r S> a, the potential V{r) can be written in the form of eq. (|57|). The string tension is 
determined by the gluon condensate and the correlation length via eq. (^). The value of 
the intercept Vq depends on the value of k and the form of the correlators D and Di . For 
K = 1 one obtains eq. (|6^). 

Therefore the linear rising of the potential at large distances as well as a negative 
intercept can be explained within the model of the stochastic vacuum. Note that the 
model of the stochastic vacuum explains only the non-perturbative part of the potential. 
The perturbative part has to be added. In the three-dimensional case it is, however, not 
easy to distinguish from the non-perturbative part because the squared gauge coupling 
has dimension of a mass. 
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